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Abstract 

This paper deals with the entanglement, as quantified by the negativity, of pure quantum states chosen at 
random from the invariant Haar measure. We show that it is a constant (0.72037) multiple of the maximum 
possible entanglement. In Line with the results based on the concentration of measure, we find evidence 
that the convergence to the final value is exponentially fast. We compare the analytically calculated mean 
and standard deviation with those calculated numerically for pure states generated via pseudorandom uni- 
tary matrices proposed by Emerson et. al. [Science, 302, 3098, (2003)]. Finally, we draw some novel 
conclusions about the geometry of quantum states based on our result. 

PACS numbers: 03.67.Mn, 02.30.Gp 



* a.datta@imperial.ac.uk 



1 



I. INTRODUCTION 



Entanglement has come to be believed as one of the cornerstones of quantum information sci- 
ence. The necessity of entanglement in quantum computation [1] and information [2, 3] tasks 
are well acknowledged. Substantial amounts of experimental effort is expended in the generation 
and manipulation of quantum entanglement. Nevertheless, the role of entanglement in quantum 
information science in general, and quantum computation in particular, is far from clear. Meyer 
has presented a version of the quantum search algorithm that requires no entanglement [4], and 
instances are known of mixed-state quantum computation where exponential speedup is attained 
in the presence of only limited amounts of entanglement [5], and other quantities have been pro- 
posed as alternate resources for the speedup [ , ]. Recent results have further illuminated the role 
of entanglement in pure-state quantum computation. It was already known, due to the Gottesman- 
Knill theorem [8], that entanglement is by no means sufficient for universal quantum computation. 
The new results [9, 10] show that, in fact, almost all pure states are too entangled to be a universal 
resource for quantum computation. Though proved in the context to measurement-based quantum 
computation, and based on the geometric measure of entanglement, which is the absolute square 
of the inner product with the closest product state, these results drive home the point that implica- 
tions on the lines of "more entanglement implies more computational power" are fallacious [10]. 
The strategy employed for proving these results can generally be termed as "concentration of mea- 
sure" [1 1], by which a typical pure state, chosen at random from the left- and right-invariant Haar 
measure, is almost always maximally entangled across any bipartition. 

Arguments based on the concentration of measure have been used to obtain average value of 
measures of correlations and entanglement in typical quantum states. Concentration of measure 
is a very powerful concept from measure theory, which puts bounds on how much the values 
of certain smooth (Lipshitz) functions can vary from their mean value. This is a consequence 
of the remarkable fact that the uniform distribution of the fc-sphere is concentrated largely 
on the equator for large k, and any polar cap smaller than a hemisphere has a relative volume 
exponentially small in k. Examples in quantum information theory include the entropy of the 
reduced density matrix, entanglement of formation, distillable common randomness [11]. The 
entropy of reduced density matrices of typical states has also been conjectured and calculated 
independently [12-15], as has been their concurrence, purity and the linear entropy [16]. Not 
much is however known of one of the most common and computable measures of entanglement, 
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the negativity [17, 18], in random Haar distributed pure states. In this paper, our endeavor will 
be to address this question. We show that the negativity of a random pure state taken from a 
Haar distribution, is a constant multiple of the maximum possible. This entanglement can also be 
generated efficiently using two qubit gates [19]. We will evaluate this constant using techniques 
similar to those in Refs. [13, 16], and confirm our results numerically using efficiently generated 
pseudorandom unitaries ["*n]. For simplicity, we will only present results for equal bipartitions, 
but the extensions to unequal splits is straightforward. 

That the negativity (defined in Eq. (2)) of random pure states is less than maximal might seem 
to contradict the statement that random pure states in large enough Hilbert spaces are close to 
being maximally entangled. This is, however, not true in general. As shown in [ i ], for a state 
residing in a Hilbert space of dimension dA x with a reduced state pa = TrB(p), and (i^ is a 
large enough multiple of log dA/e'^, then 

(l_e)-Li<p^<(l + e)-il (1) 

"A "A 

If a state satisfies Eq. (1), then its negativity is evidently near-maximal. But as the condition for 
its validity shows, this is only true when the bipartite split is quite asymmetrical. Thus, for equal 
bipartite splits, which is often of interest in quantum information science, there is no a priori 
reason to expect the negativity of random pure states to be close to maximal. This is the case 
we study here. Just to highlight the degree of asymmetry needed to have the negativity close to 
maximal, for e = 0.1 and = 2, we require ds ^ 200, and for dA = 16, ds ^ 6400. 

The outline of the paper is as follows. In Sec. II, we begin by deriving the expression of 
average negativity. It involves performing integrations over the probability simplex which are 
rewritten in terms of other nonconstrained variables, finally leaving us with a combination of 
hypergeometric functions. Sec. Ill derives the expressions for the variance in the negativity in 
terms of similar hypergeometric functions. These functions are explicitly evaluated in Sec. IV 
numerically. This is necessary as the series we have is provably not summable in closed form, 
which we discuss in brief in Appendix A. We obtain the final expression for the average negativity 
of Haar-distributed random pure states. We also compare our results with a numerical simulation 
using pseudorandom pure states generated from efficiently generated pseudorandom unitaries [20] , 
finding good agreement. We finally conclude in Sec. V with discussions about the ramifications of 
our finding on the geometry of the set of quantum states. We also discuss the prospect of extending 
the present analysis to random mixed quantum states. 
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II. NEGATIVITY OF TYPICAL PURE STATES 



The negativity is an entanglement monotone which is based on the partial transpose test of 
detecting entanglement [22] . Given a bipartite quantum state residing in Ha ® "Hb with dimensions 
jj, and called pab, the negativity is defined as 

mPAB) = MMl^, (2) 

where p^B denotes the partial transpose with respect to subsystem A, and ||cr|| denotes the trace 
norm, or sum of the absolute values of the eigenvalues of a, when a is Hermitian, as is the case 
with p^. For pure states residing in the above space, it is always possible to write a Schmidt 
decomposition [8]. This paper will only deal with the scenario p = u, the extension to the unequal 
case being tedious, but straightforward. The distribution of the Schmidt coefficients is given by 
(for p = u) [23] 

P(p)dp = iV5(l-^p,) J] {P^-P,f\{dvk. (3) 

i=l l<i<j<^^ k=l 

where is the Dirac delta function. The negativity for pure states is 
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2 Vmi (4) 



and its mean is given by 



(5) 



At the outset, it helps to change variables such that gj = rpi which removes the hurdle of integrat- 
ing over the probability simplex [13, 16], whereby 

Q{q)dq= H {q,-q,fl[e-'^^dqk = Ne~W^"-'P{p)dpdr . (6) 

The new variables qi take on values independently in the range [0, oo), and r is a scaling factor 
given by r = gj. Integrating over all the values of the new variables, we find that the normal- 
ization constant is given by = Q/T^pu), where Q = J Q{q)dq. Similarly, we find that 

J v^Q(q)dq = Q ^^py^ / VM]Pip) dp , (7) 



with r(/i) = (/i — 1)!. Notice that the first product in Eq. (6) is the square of the Van der Monde 
determinant [13, 16] 



A(q) = n - ^^•) 

l<i<j<fj, 



1 ... 1 
qi ■■■ 



ii • • • V 



Liiqi] 



Loiqf,) 
Li{qf.) 



r(/i)L^„i(gi) ... r(/i)L^„i(g^; 



(8) 



The second determinant in Eq. (8), follows from the basic property of invariance after adding 
a multiple of one row to another, and the polynomials Lk{q) judiciously chosen to be Laguerre 
polynomials [24], satisfying the orthogonality relation 



'0 

These facts in hand, we can evaluate 



j'OO 

I dqe~'^Lk{q)Li{q) = 6m . 
Jo 



(9) 



Q= f A{cO' lle-'^^dqk 

k=i 

= J2 i-'^)™Il^iTikmR{k)) dqke-'^>'LTik)-i{qk)LR^k)-iiqk 
T,_Re<S^ fc=i 



(10) 

Ri^S^ k=l k=l 

with T, R being elements of the permutation group on fx elements S^. We can now calculate the 
integral over {gi, ■ ■ ■ , g^} in Eq. (7) as 



J y^<5(q)dq 



11 fj. At 

/ VmYldqme""'- {-'^f^'^Y[^iT{m))T{R{m))LTim)^i{qm)Ln(m)-i{qn 

i,j=l ^ m=l T,ReSn m=l 



M-1 
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where 
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e-'q" Lk{q)Li{q) dq, 



(12) 



I ■ I is the determinant and we have used the orthonormality condition in Eq. (9) in the first step of 
the evaluation. We thus have 

M-l r 

(13) 
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except that the integral needs to be evaluated. 



III. VARIANCE IN THE NEGATIVITY 



Having calculated the mean of the negativity for random, Haar distributed pure states, we move 
on to calculate its variance. Based on the definition of negativity in Eq. (4), we obtain the expres- 
sion for the variance of the negativity as (cr is the standard deviation) 




(14) 



The second term has already been evaluated in the previous section, so we need con concern 
ourselves with the first term. We begin by expanding the fourth power above as 

\ 4 

A* \ A* M M M 
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2 PiPj + ^ X] Pi\/PjPk+ Y VPiPjPkPl- 



(15) 
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Each of these terms can now be individually evaluated, and omitting the details we just present the 
results as 
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IV. EVALUATING THE INTEGRALS 



Having derived formal expressions for the mean and standard deviation of the negativity of a 
random pure state, we now need to evaluate the integral in Eq. (12). To that end, we use the 
generating function for Laguerre polynomials [24] 

oo 

(1 - zY^e"'!^-^ = J2 Li{x)z^ \z\ < 1, (19) 

(=0 

and 

/" * ^ n^^l)r(...^l) ^_,_.^ f ^ . l; a . 1, 1) , (20) 

Jo n\T{a + l) V 

F being the hypergeometric function such that 

F{a, b- c-z) = }_^ — — , (21) 

n=0 ^^'"^ 

and {a)n = a{a + l)(a + 2). ..(a + n — 1) is the Pochhammer symbol. Note that if a is a negative 
integer, (a)„ = for n > \a\ and the hypergeometric series terminates. Then, 



oo «oo 

z=o 



oo 



s I e-'VLfc(a;) da; s = 1/(1 - z) 







whereby 



= s-'^T{P+l)F{^-k,f3 + l;l 

= r(/3 + 1) E E ^r^^^^7^l^(^ + (22) 

i=0 i=0 ■ 

-r(g+i)^i: (i+ffl, (23) 

i=o ^'^■^ 

and (a)„ = a(a — l)(a — 2). ..(a — n + 1) is the 'falling factorial'. Using the following identities 
for the Pochhammer symbols 

(x)„ = (-l)"(-x)„, (24) 
(-x)„ = (-l)"(x-7i + l)„, (25) 
(a;)„ = r(x + n)/r(x), (26) 



we have 




t\T{t-l + f3 + l) 



[r(t + /3 + i)F 



3F2({/3 + l,/3 + l,-A;};{l,/3 + l-/};l). 



(27) 



To get the final expression for the negativity in Eq (13), we substitute the expression for the 
integrals from Eq (27). The expressions are not very illuminating, and for the lack of an asymptotic 
expression, we present the numerical values in Table (I), and plot them in Fig. (1). See Appendix A 
for a note on the summability of the series. Anticipating a scaling in proportion to that of a 
maximally entangled state, we divide the mean expressed in Eq (13) by the maximum possible 
negativity of a /i x /i system as Mmax = (/^ — l)/2. As can be seen from Table (I), the average 
value of the negativity saturates to a constant multiple of the maximum possible. This constant is 
found numerically, and in the asymptotic limit of large n, the negativity for an equal bipartition of 
a randomly chosen Haar-distributed pure state is 



Though we have not proven this analytically, it is easily seen that the convergence is exponen- 
tial. This can be concluded from the last column in the table, which shows the difference in the 
successive values of the third column. The value of A is progressively halved as the number of 
qubits n goes up, and this shows that the negativity indeed saturates monotonically, and arguably, 
exponentially fast, to the value presented above. This is to be expected from the concentration 
of measure results [11], which means that the negativity of random states in large enough Hilbert 
spaces is close to their expectation value. 

A. Numerical verification 

As a final corroboration of our results, we test our calculations against numerically generated 
pure states. These are pseudorandom rather than random Haar-distributed. They are generated 
by applying pseudorandom unitaries presented in Ref. [20] on fiducial pure states. The negativity 
of these pure states is calculated and plotted as a histogram in Fig. (2). We compare this to an 
approximation of the cumulant generating function, and the probability distribution function for 




(28) 
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n 


A* 


{M) /Mmax 


A 


2 


2 


0.589049 




4 


4 


0.65368 


0.0646309 


6 


8 


0.686614 


0.0329346 


8 


16 


0.703378 


0.0167641 


10 


32 


0.711878 


0.0084994 


12 


64 


0.716171 


0.0042932 


14 


128 


0.718332 


0.0021611 


16 


256 


0.719417 


0.0010851 


18 


512 


0.719961 


0.0005439 


20 


1024 


0.720233 


0.0002724 


22 


2048 


0.72037 


0.0001366 



TABLE 1. Ratio of the negativity of random pure states to the maximal negativity for Haar-distributed states 
of n qubits. For an equipartition of n qubit states, fx = 2"^/^ . A is the difference between successive values 
in the third column, providing evidence for an exponential convergence of {M) /Mmax with n. 




FIG. 1. The value of the normalized negativity (TV) /Mmax of random Haar-distributed pure states. 
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FIG. 2. Distribution of the negativity of 100000 pseudorandom pure states, with n = 4 (Left) and n = 
8 (Right). The pseudorandom unitaries used were generated via the techniques of [20], with j = 40 
interactions applied for each unitary. Also plotted is the gaussian distribution function with just the first 
two moments, as given by Eq. (29), as well as the analytically calculated mean (solid vertical line) and the 
standai^d deviations (dashed vertical lines). Although the convergence of the pseudorandom construction 
of [20] to the Haar measure is not obvious, it has been shown to do so [21, 26]. 



the negativity itself P(J\f)dJ\f, given by 

P(Af)dJ\f =^ r duexp (-iMu + -^tu + dU 

J ~oo \ max J^max ^- / 

= _^e-(^-^')'/^'^"t/AA, (29) 

where M' = {M) jMmax and o' = cr/Afmax- As is evident from Fig. (2), the distributions are very 
localized, and the gaussian distribution seems quite apt. 



V. CONCLUDING DISCUSSIONS 

The negativity provides upper bounds on the teleportation capacity of a state, and its distilla- 
bility, the latter via the logarithmic negativity. It is in these two contexts that our results on the 
negativity provides new insights, not achieved through other measures. To address the telepor- 
tation capacity, the singlet distance was introduced in [18]. It is defined as closest distance any 
quantum state can get to the singlet (the ideal resource for teleportation) while undergoing only 
local operations. JVIathematically, 

A(|$),p)=inf|||$)($|-P(p)||i (30) 
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where P is the set of all local protocols, and |$) is the singlet residing in ® C"*. Note that in 
our case m = 2"/^ = 2Afmax + 1- The following result, also proved in [18] 

A(|*).,)>2(l-?^M±i) (31) 
\ m J 

then immediately leads to the conclusion that a pure quantum state I?/'), chosen at random from 
the Haar measure, will with high probability have a singlet distance given (all ~ signs here and 
henceforth apply to large n) 

A(|$), m>2(l- ^^^^l ^2(1- ^\ ^ 0.55926, (32) 



2"/2 J \ A/;. _ 

where we have used Eq. (28), which is that = 0.72037 = c < 1. This gives us a nontrivial 

-'V max 

lower bound on how close a typical pure state can be taken to a singlet by purely local operations. 
This can be recast in terms of an upper bound on the teleportation fidelity [18, 25] of random pure 
states as 

^ max ($|P(|V;)(V^|)|$) < IM±1 < ^ 0.72037. (33) 

Another application of our result can be found by using the logarithmic negativity [28] as an 
upper bound on the entanglement of distillation Ed{p)- It was shown [18] that 

Ed{p) < EAp) (34) 

where Ej^{p) = log^ \ Ip^-* | |i = log2(2A/'(p) + 1). Using this, we get (where c = 0.72037, as after 
Eq. (32)) 

(IIp^-^IIi) = c2"/2 + i_c, (35) 

whereby for a pure state \ip) chosen at random from the Haar measure, we can set the upper bound 
of distillable entanglement to be 

Eom) < log2 ((IIp"^"!!!)) ^ 2 + ^°S2C. (36) 

For the constant we present in this work, this provides us with a bound that is tighter by about half 
an ebit (log2 c ~ —0.47319). Also note that we have taken a logarithm of the average, which is 
always greater than or equal to the average of the logarithm. 

In addition to the obvious conclusions that the fidelity of teleportation and distillability of ran- 
dom pure states have nontrivial upper bounds, the above two mathematical results tell us a few 
things about the structure of the set of pure quantum states in general. Firstly, although a random 
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pure state is very likely to be highly entangled (close to maximal), it is in no way close to the sin- 
glet state, at least in trace norm. This means that a nonzero fraction of these "close to maximally 
entangled" states contain inequivalent types of entanglement which are not related by SLOCC 
operations to the canonical maximally entangled (singlet) state. A second, and probably stronger 
statement is that not only do random pure states lie in different inequivalent sets of maximally en- 
tangled states, but also that some of these classes have a greater ability to retain their entanglement 
under distillation protocols than others, resulting thereby in an overall lower distillation rate. 

This paper shows that the negativity of ra-qubit random pure states chosen from the Haar mea- 
sure is a constant multiple of maximum possible negativity, which goes as 2"/^ for an equal bi- 
partition of the state. We also provide evidence that the convergence to the asymptotic value is 
monotonic and exponentially fast. The value of the constant was not evaluated in closed form, 
and we showed why this was the case. The expression for the negativity is a sum of hyperge- 
ometric terms, and the techniques of creative telescoping show that our particular series in not 
summable. Finally, we show that the results of our analytic calculation are borne out by random 
states generated by applying pseudorandom unitaries on fiducial states. We also show that prob- 
ability distribution for the negativity is well approximated by a gaussian distribution whose mean 
and variance we obtain analytically. 

One issue that we have not addressed here is the extension of the above calculation to random 
quantum states that are mixed. This is made somewhat challenging by the fact that there does 
not exist a unique measure on the space of mixed quantum states. Since any pure state can be 
generated by applying a unitary matrix on a fiducial state, a unique measure on the space of pure 
states can be derived from that on the space of unitary matrices, which is the rotationally invariant 
Haar measure. Mixed quantum states cannot be generated in a likewise manner, and therefore, it is 
not possible to capture the distribution of mixed states via the Haar measure. However, any mixed 
state can be diagonalized by a unitary matrix, and this motivates a product measures on the space 
of mixed states Ai, which can be defined as = £ x P, where P is the usual Haar measure 
that captures the distribution of eigenvectors of the states. £ is meant to capture the distribution 
of eigenvalues, and there is no unique way of doing that. Attempts have been made [31], and 
the mean entanglement, as quantified, for instance, by the purity has been calculated, as has been 
the logarithmic negativity for tripartite states using minimal purifications [21]. The calculation of 
the negativity for states of this form will be the subject of a future publication. This will provide 
us with information about the typical entanglement(negativity) content of random mixed states, 
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which are more and more likely to be encountered as we move closer to realistic implementations 
of quantum technology. 

Appendix A: A mathematical digression 

The final expression for the negativity, though seemingly compact, is, in fact a sum of expo- 
nentially many terms. This retards the evaluation of the quantities in Table (I) drastically, unless a 
closed form is found for quantity in Eq. (27). Consequently, it would not only be interesting, but 
indeed essential to have a closed form of the above expression. For some special instances of k, I 
and (3, this is possible. Unfortunately, this is not possible for general values of k and / (this paper 
deals only with (3 = 1/2, 1). In fact, it can be shown that there exists no closed form solution for 
the sum in Eq. (27). The arguments leading to this 'tragic' conclusion are presented next. 

Theorem 1 (Zeilberger's algorithm or the method of creative telescoping [29]). Let F{n, k) be a 
proper hypergeometric term. Then F satisfies a nontrivial recurrence of the form 

J 

'^aj{n)F{n + j, k) = G{n, k + 1) - G{n, k), 

j=0 

in which G{n, k) /F{n, k) is a rational function ofn and k. 

That this theorem applies to the sum we have at hand is evident. The application of this algo- 
rithm to the expression in Eq. (27) yields third order recurrences which can be solved using the 
Gosper-Petkovsek algorithm [29, 30]. This algorithm (also called Hyper [32]) provides a com- 
plete solution to the problem in the sense that it either provides all the solution to the recurrence 
problem. On the other hand, the failure of the algorithm to come up with a solutions proves that 
the initial series cannot be summed into a closed form. It is the latter that happens in our case, 
thereby proving that the series in Eq. (27) is not summable in closed form. 

ACKNOWLEDGMENTS 

It is a pleasure to thank Colston Chandler, Anil Shaji, Adolfo del Campo and Miguel Navascues 
for several interesting discussions during the course of this work, and Martin B. Plenio for several 
comments on the manuscript. AD was supported by EPSRC (Grant No. EP/C546237/1), EPSRC 

13 



QIP-IRC and the EU Integrated Project (QAP). 



[1] R. Jozsa, N. Linden, Proc. Roy. Soc. A 459, 201 1, (2003). 

[2] L. Masanes, Phys. Rev. Lett. 96, 15050L (2006). 

[3] M. Piani, J. Watrous, Phys. Rev. Lett. 102, 250501 (2009). 

[4] D. A. Meyer, Phys. Rev. Lett. 85, 2014, (2000). 

[5] A. Datta, S. T. Flammia, C. M. Caves, Phys. Rev. A 72, 042316, (2005). 

[6] A. Datta, A. Shaji, C. M. Caves, Phys. Rev. Lett. 100, 050502, (2008). 

[7] A. Datta, S. Gharibian, Phys. Rev. A 79, 042325, (2009). 

[8] M. A. Nielsen, I. L. Chuang, Quantum Computation and Quantum Information, Cambridge Univ. 
Press, (2000). 

[9] M. J. Bremner, C. Mora, A. Winter, Phys. Rev. Lett. 102, 190502 (2009). 

[10] D. Gross, S. T. Flammia, J. Eisert, Phys. Rev. Lett. 102, 190501 (2009) . 

[11] P Hayden, D.W. Leung, A.Winter, Commun. Math. Phys. 265, 95, (2006). 

[12] D. Page, Phys. Rev. Lett. 71, 1291, (1993). 

[13] S. Sen, Phys. Rev. Lett. 77, 1, (1996). 

[14] S. K. Foong, S. Kanno, Phys. Rev. Lett. 72, 1148, (1994). 

[15] J. Sanchez-Ruiz, Phys. Rev. E 52, 5653, (1995). 

[16] A. J. Scott and C. M. Caves, J. Phys. A 36, 9553, (2003). 

[17] K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewenstein, Phys. Rev. A 58, 883, (1998). 

[18] G. Vidal, R. Werner, Phys. Rev. A 65, 032314, (2002). 

[19] R.Oliveira, O.C.O.Dahlsten and M.B.Plenio, Phys. Rev. Lett 98, 130502 (2007). 

[20] J. Emerson, Y. S. Weinstein, M. Saraceno, S. Lloyd and D. G. Cory, Science, 302, 3098, (2003). 

[21] O.C.O. Dahlsten, R. Oliveira and M.B. Plenio J. Phys. A 40, 8081, (2007). 

[22] A. Peres, Phys. Rev. Lett. 77, 1414, (1996). 

[23] S. Lloyd, H. Pagels, Ann. Phys. 188, 186, (1988). 

[24] I. S. Gradshteyn, I. M. Ryzhik, Table of Integrals, Series and Products, Academic, New York, (1980). 

[25] M. Horodecki, P Horodecki, and R. Horodecki, Phys. Rev. A 60, 1888, (1999). 

[26] J. Emerson, E. Livine and Seth Lloyd, Phys. Rev. A 72, 060302(R), (2005). 

[27] K. G. H. Vollbrecht, and M. M. Wolf, J. Math. Phys. 43, 4299, (2002). 

14 



[28] M. B. Plenio, Phys. Rev. Lett 95 090503, (2005). 

[29] M. Petkovsek, H. S. Wilf, and D. Zeilberger, A=B, AK Peters Ltd, (1997). Available for free download 

at http : / / www . math . upenn . edu/ "wilf /AeqB . html. 
[30] M. Petkovsek, J. Symb. Comp. 11, 1, (1998). 
[31] K. Zyczkowski and H-J. Sommers, J. Phys. A 34, 71 1 1, (2001). 

[32] The codes for executing this and Zielberger's algorithm are available at 
http : / / www .risc.uni-linz.ac.at/ research/ combinat / software/. 



15 



